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ABSTRACT 


The Indian subcontinent does not fall in the 
groundwave range of any LORAN-C transmission. 

As such, at present the only alternative tech- 
nique in this region for convenient and routine 
T & F intercomparison is via VLF phase measure- 
ments. At NPL, New Delhi, continuous phase 
recording of the 15 kHz transmissions from 
GBR (UK) is being made. In addition the pub- 
lished midday phase data of GBR from several 
laboratories-NPL (UK), RGO(UK) , PTB (FHG) and 
USNO (USA) - are being received regularly. In 
the present paper we discuss T & F intercom- 
parisbns between the local time scale, UTC 
(India), at NPL and those at the above mentioned 
laboratories, using the VLF phase data. A major 
factor which limits the accuracy of long term 
comparison is the seasonal variation in the VLF 
propagation delay over long paths. It is shown 
that by taking into account the seasonal propa- 
gation delay variations in a semi empirical way 
the accuracy of T & F comparisons can be 
considerably improved. In fact over a one year 
period accuracy of few parts in lO^^ in frequency 
and 1-2 fiusec in time have been obtained. 

The relative frequency offset difference between 
UTC (India) and UTC (PTB) evaluated in the present 
work as (7,0 + 0.!J)xl0“13 agrees very well with 
that obtained via the satellite experiment des- 
cribed in a companion paper. 

INTRODUCTION 

The National Physical Laboratory (NPL), New Delhi, India 

has the statutory obligation of maintaining the India 
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Standards of time and Frequency, This is being achieved 
at present with the help of- two commercial cesium clocks: 
NPL-1 (Oscilloquartz Model 3200) and NPL- 2 (Hewlett Packard 
Model 5061A with option 004), A time scale, UTC (India), 
is being maintained using the NPL-1, One of major problems 
being faced by us at NPL is that of a regular time and 
frequency transfer link between UTC (India) with those 
of the other international timekeeping laboratories. The 
Indian subcontinent does not fall in the groundwave range 
of any LORAN-C transmitter thus prohibiting the use of 
this technique. 

The first serious attempt in having an accurate link was 
made in May- June, 1979, when a two way satellite time 
transfer experiment was performed between NPL and PTB, 

West Germany (Mathur et al, 1980 and references there in). 
With the help of this experiment UTC (India) was synchroni- 
sed with UTC* (PTB), From July, 1979 onwards we are 
continuously tracking the phase of 13 kHz transmissions 
from GBR. (UK), The choice of this station was dictated by 
the fact that; (a) the received signal strength at our 
location is very good so that a precise phase tracking is 
possible and (b) there are several timekeeping laboratories 
which track this station and publish their data. We are 
regularly receiving the GBR phase data from PTB (FRG) , 

USNO (USA), NPL (UK) and RGO (UK). Utilizing the VLF data 
it has been possible to establish a fairly precise link 
between UTC (India) and UTC of the abovemen tioned labora- 
tories as will be described in the subsequent sections, 

A confirmation of the accuracy of the VLF links was made 
possible with a portable clock trip from USNO in September, 
1930, This is described briefly in a separate subsection. 

THE TECHNIQUE . 

The VLF time transfer technique has been described by 
Becker et al (1969), This basically consists of recording 
the time difference between the 1 pps of the local time 
scale and some specified phase (generally the positive 
zero crossing) of the received VIF signal appearing Just 
subsequently. This measurement is referred to as the 
•phase time', or simply 'phase' of the received signal 
relative to the local time scale. The recording is done 
at the local noon over the path mid point when VLF propa- 
gation conditions are most stable. Subtracting the daily 
values of phase measurements recorded at two laboratories 
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gives the difference between their time scales. This 
difference Is ambiguous to an additive constant which 
involves the difference between the propagation delays 
over the two paths. While this ambiguity is immaterial 
for frequency comparison between the two time scales it has 
to be eliminated if a time comparison is desired, Elimination 
of the additive constant can be achieved by an initial 
calibration v/ith the help of either a portable clock or two 
way satellite experiment. 

A ma jor factor which degrades precision and accuracy of the 
VIF time transfer is the non constancy of the propagation 
delay. There is a day to day Jitter in the propagation 
delay because of the normal variability of the ionospheric 
D-region. This is more pronounced during winter than 
during summer , as will be shown in the next section (see 
also Belrose,1953) . The Jitter evidently reduces the 
precision of the link. 'Sudden ionospheric disturbances 
(SID) due to solar flares may occasionally cause path anom- 
alies, but these generally last for short durations and 
can be isolated by careful insoection of the data (Reder, 
1971).. 

In addition, over long paths there are also systematic 
seasonal variations in the propagation delay which can 
some times be as large as 15 fi sec (lijima et al, 19S3; 
Swanson & Kugel, 1972). If these are not taken into 
account, then they may introduce significant inaccuracy in 
the time and frequency comparisons. The seasonal variation 
in the propagation delay occurs mainly due to variations 
in the D-region ionization. This consists of two parts, 
namely (a) seasonal variation of the Solar Zenith angle 
at the mid path noon and (b) variations in the mesospheric 
neutral atmospheric constituen cs, which are to sqme extent 
related to meteorological phenomena (Belrose 1963). Of the 
above, only the first part , (a), can be modeled with definite- 
ness. It has been shown by Swanson & Kugel (1972) to be 
of the form M(l-cbsX ), where, CosXis the average of the 
cosine of the solar zenith angle ^ over the path at the 
midday. M is an empirical constant dependent on the 
propagation path and to some extent on the solar activity. 

. the data 

(a) KPt (India) data : The midday phase observations of 

the 16 kHz GBR signal relative to UTC(India) recorded daily 
between 03-30 to 09-30 UT are shown in Fig. 1(a) for the 
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period 10th July, 1979 to lOth July, 1930. Twice during 
the one year time the NPL-1 underwent discontinuities in 
operation for short periods. Once during 23-25. August , 

1979 and again during 30 August - 7 Sep tember, 1980. In 
both cases it was possible to restart UTC (India) by syn- 
chronizing NPL-1 with NPL-2 and giving suitable time cor- 
rections (of 1.3 and 5.5 ysec respectively). The time 
corrections were determined from the Intercomparison data 
between these two clocks which was taken periodically 
throughout. 

We observe the following features in Fig. 1(a). The day to 
day variability is much smaller in summer than in winter. 
The rms Jitter in summer is less than 1 psec while in 
winter it is about 2.5 psec. There is a gradual drift in 
the phase data arising due to a relative frequency offset 
between UTC (India) and the transmitted frequency of GBR. 
This drift is modulated by a seasonal variation in the 
propagation delay. As mentioned in the earlier section, 
the Solar Zeni th angle related seasonal variation of the 
form M (l-c osx ) can be modelled and eliminated. To deter- 
mine cosx we follow lijima et al (1968). We divide the 
propagation path (6700km) into 10 equal segments and deter- 
mine cosX at the m idpoint of each. The average of these 
values gives c®sx over the path. To determine the value 
of M we take help of the observed diurnal variation of 
phase delay. It has been shown by Swanson and Kugel (1972) 
that around midpath noon, the t empor al variation of phase 
delay is also of the form M (1-cosX ). We thus make 
comparisons between the temporal variation of the o bserv ed 
phase delay and that of the calculated factor M (1-cosX ) 
by varying M. The value of M which gives best agreement 
between the forms of the two variations is selected. 
Utilising the whole years data the best value of M=12.5+l 
was obtained. In Figs. 2 and 3 we have shown as illustra- 
tions the observed durnal phase for gr oups of IS days in 
April and June and the calculated M(L.cosX )• The good 
agreement between the two is evident. Adopting the above 
value of M the calculated seasonal variation of M(l-cosX ) 
at the midpath hoon for the full year is shown in Fig. 1(b). 
This shows an annual variation with a summer to winter 
phase retardation of 7 psec. On subtracting out this 
seasonal variation from the observed data in Fig. 1(a) we 
get the resultant variations as shown in Fig. 4 (a). The 
phase drift now appears more regular but there still per- 
sist some short period variations during the period 
September ,1979 to February, 1980 with peak amplitudes of 
2-3 /usee. 
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(b) PTB data* The midday phase observations of GBR relative 
to UTC(PTB) is shown in Fig. 4(b). This path is only 790 km 
long and the solar zenith angle related seasonal variation 

in the phase delay is not expected to be significant. It is 
not apparent in the data also. During the period from 
September *79 to February *80, the nature of variations in the 
PTB data are very much similar to NPL (India) data in Fig. 4 
(a). It is notable in this connection that GBR-PTB path is 
almost overlapping with 1/8 of the GBR-IIPL (India) path. This 
indicates the possibility that the shorter period variations 
mentioned above are arising from this portion of the path, 

(c) RGO data* The midday phase variation of GBR relative 
to UTC (RGO) in shown in Fig,4(c). In this case also the 
path length being very short (180km) seasonal variations 
are almost absent, 

(d) USNO data * The midday phase variations of GBR Relative 
to UTC (USNO) are shown in Fig. 5 (a). The propagation path 

in this case is long (5800 km) and so the solar zenith ahgle 
related seasonal variation is significant. As we do not have 
the diurnal phase variation for this station it is not 
possible to determine M as was done earlier for NPL data. In 
absence of anything better we have just adopted the same 
value o f Msl 2.5 and calculated the variation of midday values 
of M(l-cosX), This is shown in Fig 5(b), Subtracting 
variations in Fig 5(b) from those in Fig 5(a) we get the 
resultant as shown in Fig In this case also during 

September 1979 to February 1930, we observe anomalous varia- 
tions of similar nature as In PTB and NPL (India) data. But 
these variations are of larger magnitude (-^12 ^sec peak 
to peak) 

(e) NPL (UK) data s The midday phase observation of GBR 
relative to UTC (BPL,UK) showed occasional jumps of 5,10 
and 15 psec over the whole year. These were most probably 
due to some equipment malfunction. Thus, the data from this 
station have not been used, 

RESULTS AND DISCUSSION 

For in ter comparison between UTC (India) and UTC of the 
other laboratories, the individual phase data in Figs 4(b), 
4(c) and 5(c) have to be subtracted from Fig 4(a). In 
Fig 5 we have shown t (a) UTC (India)-UTC (PTB), (b) UTC 
(Indla)r UTC(RGO) and (c) UTC (Ind ia) -UTC (USNO) . It is 
clear that of the three links the one with P-TB has the 
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minltnum jitter and anomalies due to propagation variations* 
In Fig. 6 we have drawn regression lines through the points 
v/hich represent the overall drift rates. The slope of 
these straight lines give the relative frequency offset, S, 
between UTC (India) and the other UTC scales. Following 
results are obtained. 

Sincjia,PTB = ( 7.0 + 0.1 ) x 10'^^ 

^India,RG0 = ( 7.3 + 0.2 > x 10“^^ 

^India,USNO = ( 6.9 + 0.2 ) x 10”^ 

The value of %ndia,PTB obtained above is in fairly good 
agreement with a value of (7.1 + 0.5) x 10”^ obtained via 
the satellite experiment (Mathur et al, 1980) during May- 
June, 1979, This indicates that the UTC (India) has been 
maintaining a fairly constant rate. The relative offset 
with the other two laboratories are also consistent (within 
the uncertainities) with LOHAN-C links connecting these. 
LORAN-C link results have been calculated as follows : 

SpTBjUSNO = 0.4 X 10-^ 

^PTBjRGO = 0.0 X 10-13 

To get an idea of the precision of the frequency transfer 
estimates we have made some computations of the variance, 
<r^y, on the India-PTB link which is the best of the three 
links. We get the following values : 

«^y ( X = 1 day ) = 2 X 10-H 

cr'y ( T = 10 days)= 3 x 

o"y ( c = 30 days) = 2 x 10-13 

o^y ( 100 days)= 1 x 10"^ 

The y estimate on the other two links are some what 
higher than the above. 

The graphs in Fig. 6 have not been given in absolute terms 
i.e. , the origins of the graphs are not shown. Absolute 
calibrations is necessary for actual time transfer between 
UTC (India) and other UTC scales. As mentioned earlier 
this was performed by the satellite experiment in May-June, 
1979. As a result of this experiment it was found that 
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on 29 June'79 UTC (India>-UTC(PTB) = (2,7 + 0,1) /usee. The 
VLF observations were started on 10 July onwards (as the 
GBR transmission was off for a month before this). Thus In 
?ig 6(a) we assign a value of 2,7 /asec to the ordinate 
where the regression line crosses the abscissa on 29 June' 
79, Since the difference between UTC (PTB) and UTC (RGO) 
and UTC(USNO) are. known from the LORAN-C links it is 
possible to calibrate the other two graphs also# 

PORTABLE CLOCK TRIP 

A portable clock trip from the USNO was made during 18-21, 
September, 1980, This made it possible to evaluate the 
accuracy of the time transfer accuracy of the VIF link. In 
Fig, 7 we have shown the results of only one link UTC (India) 
UTC (PTB) extended up to the end of September '80, There are 
two crosses in Fig. 7 on 29 June '79 representing the satel- 
lite experiment time transfer and on 20 September '80 
representing the portable clock result# The discrepancy 
between results of the portable clock and the regression 
line predicted by the VLF data is only 1,5 fisec which is 
quite small# 

CONCLUSION 

In this paper we have discussed a VLF time and frequency 
intercomparison link between UTC (India) and UTC of PTB, 

RGO and USNO, It is clear that the India-PTB link is the 
best of the three. It has been shown that by taking due 
account of the seasonal variations over long paths it is 
possible to. achieve frequency inter comparison to a few 
parts in 10^“^ over one year period# Time transfer can be 
achieved to an accuracy of 1-2 /isec. 

We realise, however, that this accuracy level is still not 
acceptable to BIH for inclusion of UTC (India) in the 
international UTC coordinated by them# In fact a similar 
problem must be faced by many other remote timekeeping 
laboratories which are not covered by LQRAN-C groundwaves 
or regular portable clock trips. 

Our future plans in improving our links with the inter- 
national time keeping community are t (a) reception of more 
VIJ* stations possibly CMEGA (JAPAN, Liberia, La Reunion), 

(b) having regular portable clock trips, (c) using the NNSS 
satellite signals and (d) using some geostationary satellite 
links on a regular basis# 
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1979 1980 

Fig.l (a^ Daily values of CrBTl phase relative to UTC(India). 

(b) Calculated seasonal variation factor M(l-cosX) 
for GBR-NPL path. 



Fig. 2, Diurnal variation of GBR phase for 15 days 
in A pri 1-? .3Q compared with the variation of 
M(l-cosX) (dotted line) around raidday. 
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PhoM w«iio»» et »KHi 
PrQpogotiQH Path gBR N.DaM 
101b — Z5lh jun« 1980 



ig. 3 Same as in Fig. 2 for 25 days in June' 30. 
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ig.4 (a) GBR phase relative to UTC(India) corrected 
for solar zenith angle variations, (b) G3i2 phase 
relative to U?C(PTB). (c) GBR phase relative to 
OTC (RGO). 




1979 1980 


Fig. 5 (a) GBR phase relative to irrc gT-SNO) . (b)Calculated 
seasonal variation factor M(l-cos?0 ) for GBR-US>!0 
path, (c) G3R phase relative to UGC(USNO) corrected 
*for solar zenith angle variations. 



Fig. S (a^ UTC (Indian - UTG (PTB).(b) UTC (India)- 
TITG (EGO), (c) FTC (India) - UTC (USNO). 
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Fig. 7 UTC (India)- UTC (PTB) . The two crosses 

shown- on 29 June' 79 represents the satellite 
time transfer and on 20 Sept' 30 represents 
the portable clock check. 


862 



